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. We analyze first-year data of WMAP to determine tlie significance of asym- 

Q . metry in summed power between arbitrarily defined opposite liemisplieres. We 

^ . perform tliis analysis on maps that we create ourselves from the time-ordered 

data, using software developed independently of the WMAP team. We find that 
over the multipole range / = [2,64], the significance of asymmetry is ~ 10~^, 
a value insensitive to both frequency and power spectrum. We determine the 
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O ■ smallest multipole ranges exhibiting significant asymmetry, and find twelve, in- 
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eluding / = [2,3] and [6,7], for which the significance 0. Examination of the 
twelve ranges indicates both an improbable association between the direction of 
Q ! maximum significance and the ecliptic plane (significance ~ 0.01), and that con- 

^ ! tours of least significance follow great circles inchned relative to the ecliptic at 

^! the largest scales. The great circle for I = [2,3] passes over previously reported 

2 '. preferred axes and is insensitive to frequency, while the great circle for I = [6,7] is 

to '. aligned with the ecliptic poles. We examine how changing map-making parame- 

ters, e.g., foreground masking, affects asymmetry. Only one change appreciably 
reduces asymmetry: asymmetry at large scales (/ < 7) is rendered insignificant 
^ ! if the magnitude of the WMAP dipole vector (368.11 km s~^) is increased by ~ 

I 1-3(7 (~ 2-6 km s~^). While confirmation of this result requires the recalibration 

of the time-ordered data, such a systematic change would be consistent with ob- 
servations of frequency-independent asymmetry. We conclude that the use of an 
incorrect dipole vector, in combination with a systematic or foreground process 
associated with the ecliptic, may help to explain the observed power asymmetry. 

Subject headings: cosmic microwave background — cosmology: observations — 
methods: statistical — methods: data analysis 
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1. 



Introduction 



One of the most intriguing results gleaned from the intense study of the first-year data 
of the Wilkinson Microwave Anisotropy Probe ( WMAP)"^ is the possible detection of irregu- 
larities in the temperature field of the Cosmic Microwave Background (CMB). One interpre- 
tation is that the temperature field is not a Gaussian random field, i.e., when it is expanded 
in terms of spherical harmonics 



the modes are not Gaussian random variables with zero mean and variance < Cf >, 
where 



Standard infiationary cosmology predicts a random Gaussian field. Thus non-Gaussianity, 
if it exists,^ leads to the conclusion that the standard infiationary paradigm is incomplete 
or incorrect (see, e.g., §2 of Copi, Huterer, & Starkman 2004 and references therein). How- 
ever, there are other plausible interpretations for temperature field irregularities: they may 
indicate unknown sources of foreground emission in the CMB bandpass; they may be an 
artifact of an unknown problem with the WMAP instrument; or they may be an artifact of 
the algorithm by which the temperature field maps are generated. 

A number of different approaches have been applied in the search for irregularities in the 
WMAP data. Komatsu et al. (2003) use techniques based on both the angular bispectrum 
and Minkowski functionals to determine that the WMAP data are consistent with Gaus- 
sianity, while Patanchon et al. (2004) apply a blind multi-component analysis and achieve 
the same conclusion, while at the same time finding evidence of weak residual foreground 
emission in the WMAP Q-band. However, Coles et al. (2004) examine phase correlations 
and conclude that there are departures from uniformity probably caused by the foreground, 
while Vielva et al. (2004) and McEwen et al. (2004) use wavelet-based techniques and ex- 
amine the skewness and kurtosis of wavelet coefficients, finding deviations from Gaussianity 



http : / /lambda . gsf c . nasa . gov/ product /map/m_ product s . cf m 

^See Magueijo & Mcdciros and references therein for the instructive example of non-Gaussianity detected 
in Cosmic Background Explorer data. 
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of significance^ a = 0.047 and 0.017, respectively. (The former number is McEwen et al.'s 
revised rendering of Vielva et al.'s result; Vielva et al. claim a = 0.001.) Vielva et al. exam- 
ine the frequency dependence of their non-Gaussian signal and conclude that systematic and 
foreground effects can be ruled out. Copi et al. and Land & Magueijo (2005a) use multipole 
vectors to examine the lowest multipoles (/ = [2,8]) and find correlations that are unlikely 
in Gaussian random fields at the a ~ 0.01 and 0.001 level respectively. De Oliveira-Costa 
et al. (2004) and Schwarz et al. (2004) examine the unusual alignments of the quadrupole 
and octopole planes; the latter concludes that the alignments are unlikely at the a — 10"^ 
level, and note that three of the four planes are orthogonal to the echptic plane, while the 
fourth is orthogonal to the supergalactic plane, supporting the hypothesis that unmodeled 
foregrounds are a tangible cause of irregularities in the data. Vale (2005) proposes that these 
improbable alignments arc a by-product of another foreground contamination mechanism, 
the weak lensing of the CMB dipole by local large-scale structures. 

Other researchers have looked at differences between standard hemispheres (galactic 
and ecliptic) and arbitrarily defined hemispheres. Eriksen et al. (2004a; hereafter Eriksen I) 
and Hansen, Banday, & Gorski (2004a; hereafter Hansen I) use co-added V- and W-band 
WMAP radiometer maps to find asymmetry in summed power, i.e., the ratio of summed 
powers is significantly different from one (a ~ 10~^). They find that the north pole (NP) 
direction for maximum asymmetry shifts from the galactic NP toward the galactic plane as 
the lowest I values are excluded from the analysis, and they conclude that the orientation 
of the NP at higher I values indicates that residual foreground contamination is unlikely. 
Eriksen et al. (2004b) follow up this work using Minkowski functionals and skeleton length 
to ascertain non-Gaussianity in the northern Galactic hemisphere at the a ^ 10^"* level. Park 
(2004) uses a similar mathematical framework and finds a difference between the northern 
and southern Galactic hemispheres and an asymmetry in the southern hemisphere, both 
significant at the a ~ 0.01 level. Hansen et al. (2004b) use a local-curvature method and 
find non-Gaussian features when considering the northern and southern Galactic hemispheres 
separately at the a ~ 0.01 level, and find that maximum asymmetry occurs when the NP 
is close to north echptic pole, suggesting a foreground effect. Larson & Wandelt (2004) 
examine extrema outside the most conservative WMAP foreground mask (KpO mask) and 
in addition to rejecting the Gaussian hypothesis on the whole sky, they find the variance of 
maxima and minima to be low {a — 0.01) in the Ecliptic northern hemisphere but consistent 



^Throughout this paper, significance is given as the tail integral a = PnuU, where Pnuii is a sampling 
distribution for statistic S given that the null hypothesis is true, and So is the observed statistic. Conven- 
tionally, the alternate hypothesis is accepted if a < 0.05, although the standard of acceptance is a subjective 
choice. Note that the phrases "more significant" and "less significant" refer to the value a being smaller and 
larger respectively. 
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with the null in the southern hemisphere. 

These results indicate that the evidence for irregularities in the WMAP data is tantahz- 
ing, though ambiguous, and that there is as yet no consensus as to its root cause. However, 
the hints of special alignments of hemispheres for which power asymmetry is maximized, as 
well as special alignments of the quadrupole and octopole planes, suggest that an incomplete 
knowledge of foregrounds may play a leading role. (We note the effect of one possible fore- 
ground contamination mechanism, the Sunyaev-Zel'dovich effect, was recently discounted by 
Hansen et al. 2005, in accordance with Bennett et al. 2003b.) 

One heretofore unexamined aspect of this problem is the role of the map-making algo- 
rithm itself. The authors listed above have worked with the first-year data that is provided 
by the WMAP team, and thus have not examined how the choices made by the WMAP 
team in making maps (Hinshaw et al. 2003a, hereafter Hinshaw I) affect the evidence for 
irregularities in the data. For this reason, we have independently developed map-making soft- 
ware that operates on the first-year cahbrated time-ordered WMAP data, and we analyze 
the resulting maps to determine: (a) if there is statistically significant asymmetry between 
arbitrarily defined hemispheres; and (b) if altering the map-making algorithm affects the 
observed results. 

In §2, we review the basics of map- making and discuss our own map-making algorithm, 
which we have created through the study of Hinshaw I and discussions with the WMAP tea,m. 
In §3, we analyze foreground-corrected, co- added Q-, V-, and/or W-band maps to determine 
asymmetry as a function of direction and to estimate significance both as a function of 
direction and globally over the whole sphere. In §4, we determine how the distribution 
of observed asymmetry values and the direction and significance of the maximum value 
are affected by altering the map-making algorithm. In §5, we provide a summary and 
conclusions, while in the Appendix we provide complete details of our map-making recipe. 

2. Map-Making Algorithm 
2.1. Pciradigm 

We begin by reviewing the theory underlying the making of maps for WMAP; we direct 
the interested reader to Wright, Hinshaw, & Bennett (1996), Wright (1996), Tegmark (1997), 
and Hinshaw I for more detail. The goal of map-making is to determine the maximum 
likelihood estimate of the true sky map 7^ (where the subscript p indicates data as a function 
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of sky pixel) , which is associated with the time-ordered data Tt via the relation 

ft = MTp + nt, 

where the subscript t indicates data as a function of time, M is the pixel-to-time mapping 
matrix, and fit is the vector of samples from the noise distribution. 

To simplify the problem, it is assumed that rit ~ NiQ^a^), i.e., each noise sample is 
independent and the sampling distribution is the normal distribution with time-independent 
variance cr^. (This is not actually true of the noise in raw WMAP data, and ridding the data 
of the effects of 1// noise is a major component of the data cahbration process; see §2.3.2 
of Hinshaw I.) With this assumption, the maximum-likelihood estimate is: 

fp = {M^C-'M)-\M^C-^ft) . 

One may follow Tegmark (1997) and approximate the noise covariance matrix with 
the identity matrix I, and one may furthermore assume that the matrix product M"'"M is 
diagonally dominant, so that one does not have to invert a Np x Np matrix (where Np can 
be as large as ~ 3 x 10^). (Thus in first-year data processing, the WMAP team assumes 
the beam response for each radiometer is a 5-function.) For the specific case of WMAP, the 
off-diagonal elements have magnitude 0.3% relative to the diagonal elements; this is the 
inverse of the total number of pixels that may be paired with a given pixel, given the beam 
separation of WMAP radiometers. With this assumption in place, 

where riobs.p is the number of times that sky pixel p is observed. 

2.2. Algorithm 

WMAP is comprised of ten differential radiometers covering five frequency bands: Kl 
(20-25 GHz); Kal (28-36 GHz); Ql, Q2 (35-46 GHz); VI, V2 (53-69 GHz); and Wl, W2, 
W3, W4 (82-106 GHz). The lowest frequency radiometers are meant to assist the modeUing 
of Galactic foreground emission, while the highest frequency radiometers are the more im- 
portant ones for modelling the CMB. (In this work, we concentrate upon the data of the Q, 
V, and W bands.) Each radiometer consists of two horns (denoted A and B by Hinshaw I) 
that are separated by ~ 140°. Temperature differences ATiaw.t between two points on the 
sky are measured every 1.536/n s, where n equals 12, 15, 20, and 30 for the K and Ka, Q, V, 
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and W bands respectively/ (Here, we change notation from T to AT, because our interest 
lies in anisotropics.) As described in Hinshaw I, raw temperature differences are calibrated, 
to correct for varying baselines and gains (and to correct for 1// noise). As the first-year 
data release contains only the calibrated differences ATcai,t, these are what we analyze. 

Here we summarize our map-making algorithm; more detail may be found in the Ap- 
pendix. 

We begin with a zero-temperature map and estimate ATcMB,p,true via an iterative pro- 
cess. During each iteration, the time-ordered data ATpai,* are treated sequentially. 

• For each datum, we estimate the radiometer normal vectors ua and ub- The spacecraft 
orientation as a function of time is encoded using quaternions (gi -|- q2i + q^j + q^k in 
scalar- vector form) that are recorded at the beginning of each 1.536 s science frame. 
Interpolation of quaternion elements is crucial since the spacecraft spins once every 129 
s, or fti 4° per frame. The WMAP ieam uses a Lagrange interpolating cubic polynomial 
that is anchored to the four points used in its determination. The interpolated elements 
are used to determine ua and ub in galactic coordinates, and also to determine the 
two sky pixels pa and pb that are associated with the datum. 

• The contribution of the Doppler-shifted monopole is removed from each datum: 



A7cMB,t = ATcai.t — ToX 



'P ■ {Ha - Hb) + ■ HaT - ■ nsf] , (1) 



— * 

where To — 2.725 K (Mather et al. 1999), and P is the velocity of the satellite with 
respect to the CMB rest frame expressed as a fraction of c.^ We linearly interpolate the 
barycentric velocity of WMAP, which is recorded every 30 science frames (46.08 s). We 
assume the velocity of the Sun relative to the CMB rest frame in galactic coordinates 
to be (t'|,,gab^0,gab^0,gai) = (-26.26,-243.71,274.63) km s'^ (or Vq,^,^ = 368.11 km 
s~^; based on §7.1 of Bennett et al. 2003a, hereafter Bennett I). 



''There are four data associated with each temperature difference: 

1 1 

where 3 and 4 denote hnear orthogonal polarization modes. These data are not available in the first-year 
release. (To match the notation of equation 1 of Hinshaw I, replace A and B with 1 and 2, respectively.) 

^Thc effect of the second-order term, which the WMAP team does not include in their map calculation, 
is < 2/xK. See Figure 3. 
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Running sums are continuously updated to estimate A7cmb,p during the {i + 
iteration: 

Arp _ ATcMB.t + (1 - Xirn)^TcMB,pB,i /^n 

^-'CMB,pA,i+l — ■^-i CMB,pA,j+l + 7~ \'^) 

i- 1 Xiia 



and 



_ I ... -^TcMB,t + (1 + a;im)ATcMB,pA,i /rN 

^-l-CMB,pB,i+l — ■^-'CMB,ps,i+l + , [O) 

1 Xii 



where loss-imbalance parameter tabulated for each radiometer horn (~ 10 

10~^; Table 3 of Jarosik et al. 2003), and Wt is a statistical weight based on instrument 
temperature that is expected to differ from 1.0 by only ± 1% (G. Hinshaw 2004, private 
communication). We assume Wt — 1. 

There are three situations in which the running sums are not updated. First, if the 
data are flagged as bad, or second, if there is a planet within ^cut degrees of either 
radiometer A or B, then nez^/ier running sum is updated. (The VFMAPteam assumes 
^cut — 1-5°.) Third, if either radiometer is pointing to within the Galactic mask the 
running sum for the opposite radiometer is not updated. This is to mitigate the effect 
of elliptical radiometer beams, which cause the differential signal associated with bright 
sources to change with spacecraft orientation and thus to bias map estimates (Hinshaw 
I, §3.3.4). 



After all the time-ordered data are treated, the weighted average of ATcMB,p,i+i is 
determined by dividing by nobs,p; this estimate is passed to the next iteration. 

The iterative process determines the spatial distribution of relative temperature dif- 
ferences, and does not determine a zero point. As noted in the Appendix, the WMAP 
team determines a zero point via model fltting. We assume the WMAP zero-point values, 
i.e., we apply a constant offset to each pixel such that our mean temperature matches that 
determined by WMAP. 



2.3. Example Maps 

We display combined Q-, V-, and W-band radiometer maps in Figure 1. Maps are 
flrst created for each individual radiometer, then combined using the weighting scheme of 
equation (24) of Hinshaw et al. (2003b; hereafter Hinshaw 11) . 



-8- 



While our maps compare very favorably with their counterparts shown in Figure 2 of 
Bennett I, they do not match exactly. In Figure 2, we show the differences between combined 
Q-, V-, and W-band WMAP maps and our combined maps. Two residual structures are 
apparent in this figure. The first is a speckled pattern aligned with the scan pattern for the 
first day of observations. Its cause is unknown: while it may be a thermal effect manifested 
as significant departures (~ fO%) oiwt from unity, a timing mismatch between datascts has 
also been suggested (Hinshaw, private communication). The fact that the values of Wt are 
not included in the data release makes it more difficult to determine the cause; we would 
ask that these values be included in future data releases. In this particular instance, the 
effect of the missing data is not grievous; we have examined difference maps made with and 
without the first day's data and determined that the speckled pattern has negligible effect 
on maps and asymmetry values. Since our bias is to avoid arbitrary data cuts when possible, 
we include the first day's data in all analyses discussed below. 

The second residual structure is a dipole that has amplitiide 7//K and direction 
(^gaij^'gai) ~ (—86°,— 19°) in all three bands, as determined by the HEALPix^ IDL routine 
remove_dipole.^° The appearance of residual dipole structures is not surprising. When 
calibrating the data, the WMAP team assumed the barycentric velocity vector estimated 
from CO^EDMRdata, (^^0,gai, ^^0,gai> ^^agai) = (-24.57,-244.47,275.00) kms"^ (G. Hinshaw, 
private communication; based on Bennett et al. 1996, in which the vector elements are not 
explicitly stated) , and in Bennett I it is mentioned that residual dipole structures are observed 
and are used to improve estimates of the dipole temperature and direction. We conclude 
that the residual dipole that we observe reflects a combination of numerical differences the 
dipole unit vectors and statistical weights that we use and those assumed by the WMAP 
team. For completeness, we determine that removing the speckled pattern by not using the 
first day's data has negligible effect on the residual dipole amplitude (< 0.5%). 

In Figure 3, we show the difference maps of Figure 2 with the residual dipoles removed. 
Another structure, aligned with the CMB dipole, is apparent in these figures; this is the effect 
of including the second-order term when computing the Doppler shift of the monopole, which 
the WMAP team does not include in their calculations. 

Histograms of the temperature differences in all pixels are displayed in Figure 4. If we 
examine only those pixels for which riobs differs by more than 1% between our maps and 



^The HEALPix software is described in Gorski, Hivon, & Wandelt (1999) and is available at 

http ; / /www . eso . org/science/healpix. 

^°We specify galactic longitude I and latitude h with a subscript "gal" so as not to confuse longitude with 
the multipole number I. 
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those of WMAP (presumed to be "speckled pixels"), we find that in the Q band, the fraction 
of pixels which are speckled pixels at 10, 20, and 40 fiK is ~ 0.07, 0.25, and 0.70 respectively. 
In the V band the respective fractions are 0.08, 0.20, and 0.70, while in the W band, they 
are 0.03, 0.10, and 0.40. Speckled pixels thus dominate the tails of the histograms. 



3. Determination of Asymmetry in Generated Maps 

To determine asymmetry in summed power between arbitrarily defined opposite hemi- 
spheres, we follow Eriksen I and Hansen I by principally analyzing a combined map consisting 
of foreground-corrected V- and W-band maps. We perform foreground corrections using the 
method of Bennett et al. (2003b; hereafter Bennett 11), and combine the maps using equa- 
tion (24) of Hinshaw 11. To compute cut-sky pseudo-C/ power spectra, we implement the 
Master algorithm of Hivon et al. (2002), summarized in Appendix A of Hinshaw 11.^^ This 
differs from Eriksen I and Hansen I, who use the maximum likelihood method of Hansen 
et al. (2002) and Hansen & Gorski (2003) to compute cut-sky power spectra. As discussed 
in §2.2 of Efstathiou (2004), use of the Hivon et al. method is not optimal in that there is 
a estimator-induced variance that makes it an inaccurate determiner of power in the \ow-l 
regime, with the variance of the sampling distribution increasing with mask size. In this par- 
ticular work, however, the use of the Hivon et al. method does not have a detrimental effect 
in that we are not interested in power amplitudes or the magnitudes of the ratio per se, but 
rather in how the observed ratios compare with sampling distributions that are determined 
via simulations and thus have the estimator-induced variance built in. 

To compute asymmetry as a function of location on the sky, we define 101 indepen- 
dent NP orientations on the northern galactic hemisphere (/gai = [0°,324°] with A/gai = 36° 
and 6gai = [0°,90°] with A6gai = 9°; /gai is not defined for 6gai = 90.0°). We apply two 
Kp2-|- hemispherical masks for each NP orientation and compute power spectra for the north- 
ern and southern hemispheres for Z < 64. We then calculate asymmetry over arbitrary 
multipole ranges, using arbitrary weightings (e.g., wi — 1, or wi — l{l -\- 1)), via: 

^obs(^gab ^gab ^miiu ^max) '^l) — ; „ott • (6) 



^^We compute spherical harmonic transforms using the ccSHT hbrary; see 

http : //crd. Ibl . gov/~cmc/ccSHTllb/doc/. 

^^Our grid is not an equal- area grid on the sky (unlike Eriksen I and Hansen I) and is primarily designed 
to assist visualization via contour plots. 
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Fig. 1. — From top to bottom, combined Q-, V-, and W-band radiometer maps created by 
making maps for each individual radiometer and summing them using the weighting scheme 
of equation (24) of Hinshaw II. 
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Fig. 2. — From top to bottom, combined Q-, V-, and W-band radiometer difference maps 
(ATwMAP — AT) created with combined WMAP maps and the maps shown in Figure 1. 
There are two structures present in these data: a speckled pattern related to the statistical 
weighting of first day data; and a dipolar structure of amplitude ~ 7 /xK in the direction 
(/gai,&gai) ~ (-86°, -19°). See also Figure 4. 
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Fig. 3. — Same as Figure 2, but with the residual dipolar structure removed. The remaining 
structure is due to our including the second-order term when computing the Doppler shift 
of the monopole; the WMAP team does not include this term in their calculations. See also 
Figure 4. 
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Difference Histogram (Q-band) 



Difference Histogram (V-band) 





AT {fjK) 



AT {fjK) 



Difference Histogram (W-band) 





AT (/iK) 



Fig. 4. — Difference histograms ATwmap — AT for the Q-, V-, and W-band data. The sohd 
hne is a histogram of the data shown in Figure 2, while the dashed hne corresponds to Figure 
3. 
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(In this work, we assume wi — l{l + 1).) We examine each direction separately because of 
the number of pixels within the north and south masks is generally unequal; this causes the 
expected mean asymmetry value derived from simulations to differ from unity and to vary 
as a function of direction. 

We compare our results with simulated samphng distributions. We use the HEALPix 
facility synf ast to generate random universes by sampling from a given power spectrum. 
Current theoretical bias leads us to adopt the LCDM power spectrum derived from a running- 
index primordial spectrum using WMAP, CBI, and ACBAR data (which we denote RI). We 
simulate 512 random universes, generating two copies of each. We convolve the two copies 
with the 21' V-band and 13' W-band beam responses respectively. From the two band- 
specific maps, we generate six noisy radiometer-specific maps, where we assume Gaussian 
pixel noise ~ N{0,ao/^ynobs^) a-nd use tabulated radiometer-specific values of (Jq.^^ We 
co-add and analyze the six noisy simulation maps in the manner of our six observed maps, 
generating a distribution of 512 asymmetry values in each defined direction (Zgai, &gai)- 

We observe that the sampling distribution for ^(/gai, &gai); the ratio of two Gaussian- 
distributed random variables under the null hypothesis, is approximately lognormal: 



ln(A-/x)^ ^' 

\/2(T 



(7) 



This is expected because if, for instance, the distribution mean is unity, one would expect 
the same probability for observing asymmetry values of, e.g., 0.5 and 2.0. We fit lognormal 
distributions to the data to estimate significances (see Figure 5); otherwise, estimates would 
be limited to values a > 1/512 ~ 10^'^ and would be inaccurate in the distribution tails. 
For a given set of parameters (/gai, 6gai, ^min, ^max), we sum over the multipole range, obtain 
best-fit values fio and cXo by maximizing the sum of the log-likelihoods for each datum, and 
determine what we call the directional significance, aAix'- 

Qidir = min[ / dAf{A,iio,cro) , / dAf{A,iio,cro)]. (8) 
Jo JAob^ 

ctdir is an estimate of the probability that a random Gaussian field exhibits at least as much 
power asymmetry as observed in the data in the specific direction (Zgai, &gai)- 

While the directional significance is useful for, e.g., visualization, it is not sufficient to 
establish the global significance of asymmetry for a given dataset, agiobai, which is an estimate 
of the probability that the random Gaussian field exhibits at least as much asymmetry as 



13 

http : / /lambda . gsf c . nasa.gov/product/map/pub_papers/firstyear/supplement/WMAP_supplement . pdf 



- 15 - 



observed in the data in any direction. We cannot easily obtain ctgiobai because the directional 
distributions are correlated, with the degree of correlation dependent on the multipole range 
considered. (We observe that because the distributions are correlated, the simplest Bonfer- 
roni correction— multiplying the smallest observed value of directional significance, adifmim 
by the number of examined directions, 101— is too conservative.) A rigorous approach would 
involve transforming our simulation results to log-space and numerically determining the 
correlations between all combinations of grid points, which would allow us to specify the 
101-dimensional Gaussian probability function. We would then examine this function using 
methods described in, e.g., Adler (2000) to estimate ctgiobai- However, as this calculation is 
very computationally intensive, we adopt a simpler approach: we compare log(adk%in) "with 
its probability distribution, which we find to be approximately an extreme- value distribution: 



H - l0g(Q;dir,min) f 1^ - log(Q;dir,mmy 

^ ■ — - — exp ' 



l0g(adl'r%in) ~ 5[log(adir,min)] = ^ exp 

' ^ (9) 

As this distribution is the limiting distribution for the largest element of a set of independent 
samples from a normal distribution, this finding is not unexpected. In a similar manner to 
that described above, we determine best- fit values /ig and Uo (see Figure 5); the estimated 
global significance is then: 

Vo - log(a2lf 



ttgiobai = 1 - exp 



exp 



(10) 



We consider Ogiobai < 0.05 to be significant. We make no claim that this computation is the 
optimal proxy for a full Gaussian random field calculation; more accurate and/or powerful, 
yet still computationally simple, hypothesis tests may exist. 

In Figure 6 we show the direction of maximum significance and CKgiobai for the multipole 
range I — [2,iinax]- Unlike Hansen I, who find that the direction of maximum significance 
is always toward the galactic NP if Zmin = 2, we find that it moves from the vicinity of 
the galactic NP toward the galactic plane as /max increases. For instance, for the range 
/ = [2,40], we observe that a^!^"'' ~ 10~^Q;dir • For /min > 5, we find results consistent with 
those of Hansen I. Over the full multipole range that we consider, / = [2,64], we find a 
maximum directional significance a^ir — 3.3x10"^ in the direction (/g^i, 6gai) = (72°, 9°), with 
global significance ccgiobai = 2.6x10"^. (This is consistent with the ~ 10~^ value determined 
by Hansen I using 2048 simulations, without fits of samphng distributions.) This result 
indicates clearly that there is power asymmetry in the WMAP data. 

At what scales is power asymmetry most evident? Mindful that there could be a number 
of causes manifesting themselves at different scales, we determine the smallest ranges with 
significant asymmetry. We compute ctgiobai in multipole bins with minimum width Al — 2. 
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1 2 3 4 5 



Fig. 5. — Example fit of a lognormal distribution to 512 simulated asymmetry values for 
given direction (/gai, &gai)) which is used to estimate aair (top panel; see equation 8), and 
example fit of an extreme- value distribution to 512 simulated values of adir.min, which is used 
to estimate agiobai (bottom panel; see equation 10). 



-17- 




max 



Fig. 6. — Direction of maximum asymmetry in galactic coordinates (top two panels), and 
the global significance agiobai (bottom panel), computed in bins I = [2, /max] for the RI input 
power spectrum. The horizontal dotted line in the top two panels indicates the position of 
the ecliptic NP, while the dot-dashed line in the bottom panel marks ctgiobai = 0.05. 
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(We do not examine each multipole singly because the computation of ctgiobai is comphcated 
at low / by the fact that the sampling distributions for C; are truncated at zero, resulting in a 
considerable number of asymmetry values that are either zero or infinite. This is not an issue 
for A/ > 2.) We start with Al = 2, and increase the bin size until no more significant bins are 
found. For a given Al, we examine all bins (e.g., for Al = 2, we examine / = [2,3], then [3,4], 
etc.), and mark those for which agiobai < 0.05. If two or more significant bins overlap, we 
mark only the most significant bin. For each Al, we examine only those multipoles that have 
not yet been marked as part of a significant bin (e.g., if we find that only the bin I — [2,3] is 
significant for Al — 2, the first bin we would examine for AZ = 3 is [4,6], and not [2,4]). We 
note that this binning algorithm is not meant as a rigorous statistical exercise (we examine so 
many bin combinations that a number of of the chosen bins may be false positives) but rather 
as a guide to further examination of the data. In Figure 7 we show the result: we find twelve 
significant ranges, including two pairs and one triplet of (nearly) contiguous bins (separated 
by at most one bin). The gaps between significant bins arc treated as single, insignificant 
bins for purposes of plotting. As shown in the bottom panel of Figure 7, combining (nearly) 
contiguous bins does not necessarily yield a more significant result. The large number of 
significant bins and the fact that we discern no general trend in the direction of lead 
us to conclude that the process(es) contributing to the observed asymmetry (e.g., residual 
foreground emission) are spatially complex. 

The bins I = [2,3] and I — [6,7] appear anomalously significant, with ctgiobai ^ 10~^ (the 
value depends very sensitively upon the fit parameters of the extreme- value distribution and 
we cannot determine it accurately). For / = [2,3], the direction of maximum significance is 
(^gai,&gai) = (144°,9°), with C2^3 = {9.7,1.1} and C^g = {301.7,755.9}, while for / = [6,7], the 
corresponding values are (/g^i, fegai) = (180°,27°), C^^j = {303.6,196.7}. and C^j = {51.7,0.0}. 
We confirm that the simulated values of ctgiobai in each case follow extreme- value distributions 
with no discernible deviation that may affect the estimation of significance. Thus we take 
the observed significances at face value and conclude that there is highly significant power 
asymmetry at the largest scales. We note that we determine the direction of maximum 
asymmetry for I — [2,3] to be in the galactic plane, while for Hansen I it is the galactic NP 
(for I — [2,4]; we confirm that our result does not change if we include I — A). 

To determine the sensitivity of ctgiobai to changes in frequency, we simulate separate sets 
of co-added Q-, V-, and W-band data, as well as a set of co-added data from all three bands. 
In all cases, we find results that are similar, both qualitatively and quantitatively, to those 
presented above; for instance, over the multipole range I — [2,64], we find ccgiobai — 3.0x10"^ 
(Q), 6.2x10"^ (V), 7.1x10"^ (W), and 2.3x10"^ (all bands). This finding is consistent with 
that of Hansen I, who analyze Q-band data with a KpO mask and find persistent asymmetry. 
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We also repeat our analysis using four alternate theoretical^ and phenomenological 
power spectra as input to simulations: 

1. a running- index spectrum similar to the RI spectrum, with 2dF and Lyman-a data 
also included in the computation (BF); 

2. the theoretical LCDM power spectrum derived from a power-law primordial spectrum 

(PL); 

3. the observed WMAP "One- Year Combined TT Power Spectrum" (TT);^^ 

4. and the same as immediately above, but optimally smoothed via local linear smoothing, 
with bandwidth chosen by leave-one-out cross-validation (see, e.g.. Fan & Gijbels 1996; 
OP). 

For each spectrum, we compute ttgiobai within the previously defined RI bins, as opposed to 
defining a new set of bins. This has little effect on our results, which we display in Figure 
8. We find that ctgiobai is largely insensitive to the choice of power spectrum; while the 
variance of the distribution of simulated asymmetry values differs for each, the differences 
are insufficient to qualitatively affect results. However, we do note that the range I — [49,51], 
significant given the RI spectrum, is not significant for any other spectrum. 

In Figures 9a-b, we display contours of ctdir for each of the twelve globally significant 
multipole ranges shown in Figure 7. The direction and significance of maximum asymmetry 
in each panel in these figures maps directly back to the values plotted in the top three 
panels of Figure 7. Interestingly, these contours exhibit evidence of a relationship between 
the ecliptic plane and both the direction of maximum significance and the contours of least 
significance (ctdir ~ 0.2-0.5). We find that in four of the twelve plots, the direction of 
maximum asymmetry lies within ^ 4° of the ecliptic plane and in eight of twelve, within 
20°. The coarseness of our coordinate grid precludes an exact quantitative interpretation, 
but we note that the probability of sampling a point within 4° of the ecliptic plane is 0.07, 
and within 20°, 0.34; the probability of what we observe is Pbinomiai i$ 0.01 in both cases. 

We also observe that, for large scales, the contours of least significance appear to follow 
great circles; furthermore, the great circle for I — [6,7] passes over the echptic poles. (Great 
circle-like structures appear at other scales as well, but not as strikingly.) While the signifi- 
cance of such an alignment of contours is difficult to quantify, we visually examine contour 



http : / /lambda . gsf c . nasa.gov/product/map/wmap_lcdm_models . cf m 
'http://lambda.gsfc.nasa.gov/data/map/powspec/map_comb_tt_powspec_yrl_vlpl.txt 



-20- 



plots created from simulated datasets and find no discernible trend toward the formation 
of great circles. In Figure 10, wc concentrate on the lowest four multipole ranges shown in 
Figure 8. The dotted and dash-dotted lines in these figures represent great circles inclined 
relative to the ecliptic and galactic planes, respectively. Note that we place the lines by eye 
and make a purely qualitative comparison of goodness-of-fit. For / = [2,3] and / = [6,7], great 
circles defined in ecliptic coordinates with inclinations ~ 30° and ~ 90° respectively appear 
to fit the contours of least significance better than their galactic coordinate counterparts; for 
I — [8,10] and [11,14], there is no clear preference for either coordinate system. In Figure 11, 
we demonstrate that the contours and great circles for I — [2,3] are frequency independent: 
they are observed within each radiometer band separately, with the results in each band 
relying on completely independent sets of simulations. We conclude that whatever causes 
power asymmetry at the largest scales may be related to the ecliptic, but at the same time 
the underlying process is frequency independent. 

In Figure 12, we show the relationship between the great circle for the / = [2,3] and the 
axes defined by de Oliveira-Costa et al. that maximize the angular momentum dispersion of 
the quadrupole and octopole in their wave-function paradigm, as well as the four normals to 
the planes defined by the multipole vectors of Schwarz et al. All except the w*^^'^'^) vector 
of Schwarz et al. point in the same general direction (the "axis of evil," as dubbed by Land 
& Magueijo 2005a), near a contour of minimum significance (ctdir > 0.4). However, they 
also point near the CMB dipole. In the next section, we examine the effect of changing the 
magnitude of the velocity of the CMB dipole on observed asymmetry. 

4. Effect of Map-Making Algorithmic Assumptions Upon Observed 

Asymmetry 

A heretofore unexamined aspect of the asymmetry issue is the effect of altering the 
map-making algorithm. For instance, what if maps were made using no foreground mask at 
all, or a reduced mask, rather than the Kp8 mask with no edge smoothing? Because of the 
elliptical beam response, maps would be affected by propagated echoes of bright Galactic 
sources, as noted by Hinshaw I. But would these echoes have any discernible effect on the 
observed asymmetry? (We realize that in this example using no foreground mask at all is 
clearly a wrong step no one would make; our goal in these tests is simply to establish how 
sensitive asymmetry is to changes in the map-making algorithm.) 

The alterable aspects of our map-making algorithm, both major and minor, fall into a 
number of categories: 
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Fig. 7. — Direction of maximum asymmetry in galactic coordinates (top two panels), the cor- 
responding significance a^ir (third panel) , and the global significance ctgiobai (bottom panel) , 
as a function of /, for the RI input power spectrum. See §3 for a description of how bin 
widths are determined. The horizontal dotted line in the top two panels indicates the posi- 
tion of the ecliptic NP, while the dot-dashed line in the bottom panel marks Ogiobai = 0.05. 
The additional dotted lines in the bottom panel indicate the values of ctgiobai if we combine 
(nearly) contiguous significant ranges: / = [6,10], [6,14], [8,14], [22,25], and [43,51]. 
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Fig. 8. — Similar to the bottom panel of Figure 7, showing global significance agiobai as a 
function of / for the five different input power spectra listed in the text. The top panel of 
this figure is the same as the bottom panel of Figure 7. This figure demonstrates that CKgiobai 
is relatively insensitive to the choice of power spectrum. 
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Fig. 9a. — Directional significance contour plots showing the significance of observed asym- 
metry as a function of NP location in (Zgai, 6gai) within multipole bins in which the global 
significance of asymmetry is a < 0.05. The minimum bin width is AZ = 2. The direction of 
maximum significance is labeled with a triangle, the cross and solid line indicate the ecliptic 
NP and plane respectively, and the asterisk and the dashed line indicate the supergalactic 
NP and plane respectively. 




Fig. 9b. 
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Fig. 10. — Directional significance contour plots for the lowest four multipole bins shown 
in Figure 9a (see this figure for an explanation of symbols). All four plots exhibit contours 
of least significance that appear to follow great circles. The dotted and dash-dotted fines 
represent great circles inclined relative to the ecliptic and galactic planes, respectively. The 
lines are illustrative and their placement is done by eye. 
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Fig. 11. — Same as the upper left panel of Figure 10, with co-added Q-, V-, and W-band data 
displayed in the upper right, lower left, and lower right panels respectively. See Figures 9a 
and 10 for explanations of symbols. We conclude that the process causing power asymmetry 
at the largest scales is insensitive to frequency. 
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Fig. 12. — Upper left panel of Figure 10 with the great circle in galactic coordinates removed. 
The diamonds represent the pointing direction of the quadrupole (lower right) and octopole 
(upper left) preferred axes defined by de Oliveira-Costa et al., while the squares represent the 
pointing direction of normals to the planes defined by multipole vectors derived by Schwarz 
et al. The filled circle shows the direction of the CMB dipole. 
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1. The mapping paradigm, involving issues which are too complex to address in the 
current work: 

a. The assumption of Gaussian pixel noise and its subsequent use in least-squares 
estimates. 

b. The assumption that the matrix product (M^M)^^ is diagonally dominant and 
is ~ n~^g. This goes hand- in-hand with the assumption that the radiometers have 
(5-function spatial response, an assumption which we expect has little to no effect 
in the low-Z regime. 

2. Aspects which are currently not alterable given the data we have: 

a. The algorithm by which the raw data are cahbrated (Hinshaw I, §2.2). 

b. The assumption of equal statistical weights for each datum. 

3. Aspects which we directly examine: 

a. The assumption that To — 2.725 K. 

b. The assumption that the Sun's velocity relative to the CMB is (^©^gai) '^©.gai' '^0,gai) — 
(-26.26,-243.71,274.63) km s-i (or ^;0,gai = 368.11 km s-^). 

c. The effect of including the loss- imbalance parameters (Table 3 of Jarosik et al.). 

d. The use of the Kp8 galactic mask without edge smoothing. 

e. The use of Lagrange interpolating polynomials to determine radiometer normal 
vectors as a function of time. 

f. The assumption that the planetary cut is ^cut = 1-5°. 

g. A step beyond map-making: the use of the same foreground maps, in the same 
proportion, as the WMAP team (Bennett II, §6). 

We also examine the sensitivity of asymmetry to temperature map zero points, which the 
WMAP team derives by fitting a phenomenological model to southern galactic hemisphere 
data (Bennett II, §5). Normally, changing the zero point, i.e., changing the monopole com- 
ponent of the anisotropy map, should have no effect on higher multipoles. However, power 
leakage from the monopole in, e.g., the Master algorithm means that any uncertainty in the 
zero point can in theory affect measurements of asymmetry. The average temperatures in 
the V-band and W-band maps are 78 /uK and 75 /xK, respectively, with variance ~ 1 /xK. 
We confirm that changing the zero point by amounts up to 10 fiK has negligible impact upon 

•^global • 
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After altering the value of a given parameter (while holding the other parameter values 
to their default values), we create new V- and W-band radiometer maps. We combine 
these maps into a new co-added map that we analyze using the methods described in §3. 
To determine if altering mapping parameters can reduce or eliminate all asymmetry, we 
concentrate upon the range / = [2,64] (Table 1); however, we also examine the range / = [2,3] 
to gauge the effect of changes on large-scale asymmetry. 

Monopole Temperature. We find that altering the monopole temperature by la (0.001 
K) has little effect on the distribution of significances, and none on ctgiobai- 

Dipole Velocity Magnitude. We test changing the magnitude of the dipole velocity 
along the direction (/gai,fegai) = (263. 85°, 48. 25°). (Testing changes on a grid of dipole locations 
is currently too computationally intensive.) We find that positive changes ^ lo" (> 2 km 
s~^) can have a marked effect upon ctgiobai at the largest scales, even sometimes rendering 
asymmetry insignificant (Figure 13). At smaller scales, the effect on asymmetry is minimal 
(e.g., bottom panels of Figure 13). Note that the velocities that we assume in map-making 
do not match that used to calibrate the data (the COBE DMR dipole), so we cannot state 
with certainty that changing the dipole velocity can by itself eliminate asymmetry at the 
largest scales. 

To determine the effect of using the Master algorithm upon quadrupole power, we ex- 
amine both uncorrected and corrected cut-sky power spectrum estimates. As illustrated 
in Figure 14, the uncorrected Kp2+hemisphere cut-sky power spectrum estimates for the 
quadrupole (triangles and boxes for the north and south hemispheres respectively in the 
middle panel) are greatly affected when we increase the dipole velocity, with the ratio of 
estimates tending toward unity (thus reducing the significance of asymmetry). The subse- 
quent shifting of dipole power to the quadrupole via mode coupling in the Master algorithm 
changes the quadrupole power by < 10% (varying slightly as a function of direction). We 
conclude that the observed reduction in the significance of asymmetry is primarily due to 
increased quadrupole power in each hemisphere, and is not purely an artifact of the Master 
algorithm. 

It is important to note that increasing the sky coverage reduces the effect of dipole 
velocity upon quadrupole power, and if we use the Kp2 mask alone, the effect is minimal 
(circles, middle panel). Thus changing the dipole velocity has neghgible impact on "full sky" 
power spectra, and does not help to explain any supposed irregularities in these spectra (see, 
e.g.. Land & Magueijo 2005b and Vale 2005). It is also important to note that the increase 
in the dipole power for the Kp2 cut sky at > Icr (top panel. Figure 14) does not necessarily 



Table 1. Effect of Changing Mapping Parameters on Power Asymmetry for 1 = [2,64] 



Changed Parameter 


ccdir Quartiles 


'3Jdir,min 


(1 h \ 

I, 'gal, 'J gal) 


global 




Default Map 


(0.152,0.0132,2.56x10-3) 


3.30x10-6 


(72°,9°) 


2.55x10" 


-4 


No First Day Data 


(0.152. 0. 0132.2. 58xl0--'^) 


3.30x10-6 


(72°,9°) 


2.55x10" 


-4 


To = 2.724 K (-Iff) 


(0.147,0.013,2.58x10^3) 


3.30x10-6 


(72°,9°) 


2.55x10" 


-4 


To = 2.726 K (+1ct) 


(0.152,0.013,2.58x10-3) 


3.30x10-6 


(72°,9°) 


2.55x10" 


-4 


Vd = 366.24 km s"^ {-la) 


(0.144,0.013,2.26x10-3) 


4.59x10-6 


(72°,9°) 


3.46x10" 


-4 


Vd = 369.98 km s"^ (+1(7) 


(0.147,0.014,2.75x10-3) 


2.35x10-6 


(72°,9°) 


1.88x10" 


-4 


Vd = 371.85 km s"^ (+2(7) 


(0.157,0.017,3.77x10-3) 


2.10x10-6 


(72°,9°) 


1.69x10" 


-4 


Vd = 373.72 km s"^ (+3(7) 


(0.147,0.024,5.07x10-3) 


2.63x10-6 


(72°,9°) 


2.08x10" 


-4 


No Mask 


(0.148,0.013,2.86x10-3) 


4.11x10-6 


(72°,9°) 


3.13x10" 


-4 


No Interpolation 


(0.093,7.20 X 10-3,4.39x 10-^) 


5.70x10-8 


(72°,9°) 


6.20x10" 


-6 


Linear Interpolation 


(0.152,0.0132,2.58x10-3) 


3.30x10-6 


(72°,9°) 


2.55x10" 


-4 


Planet Cut: ^cut = 0.0° 


(1.45X 10-^5.77x 10-l^2.17x 10-^9) 


1.40x10-26 


(72°,9°) 


0.0 




Planet Cut: ^cut = 0.5° 


(0.147,0.0154,2.28x10-3) 


2.63x10-6 


(72°,9°) 


2.08x10" 


-4 
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rule out higher velocities, as the offsets of the tested velocity from the true velocity and from 
the COBE DMR velocity will have a non-negligible effect on the dipole power estimate, and 
the contribution of the primordial dipole component is unknown but could be significant (~ 
1000 ijK). 

In Figure 15, we display how the contours of directional significance evolve as the mag- 
nitude of the dipole velocity increases. In addition to an increase in Ogiobab we observe that 
the direction of maximum significance migrates toward the ecliptic plane, with the great 
circle still evident, though not as strongly, at 3cr. 

Loss Imbalance. We set in eqs. (3) and (5) to zero. This has no discernible effect on 
asymmetry. 

Masking. To gauge the effect of masking, we make maps with no mask at all. This has 
virtually no effect on asymmetry, presumably because the scale size of the aforementioned 
galactic echoes is smaller than those probed in our analysis. 

Interpolation to Determine Spacecraft Orientation. The quaternions that encode 
spacecraft orientation arc recorded at the beginning of each 1.536 s science frame, while 
data are sampled up to 30 times per frame. Interpolation is normally done using Lagrange 
interpolating polynomials (see the Appendix for details). We test the effect of making 
maps with no interpolation and with linear interpolation. Turning off interpolation leads to 
pointing errors on the order of degrees, and asymmetry becomes more pronounced, with the 
significance increasing by two orders of magnitude (ctgiobai ~ 10~^ for I = [2,64]). The use 
of linear interpolation has virtually no effect on asymmetry. This indicates that the WMAP 
team's use of Lagrange interpolating polynomials is a conservative algorithmic choice, and 
that the small errors in the recorded orientation of the spacecraft and the determination of 
instrumental boresights (< 1'; see §3.3.1 and 3.3.2 of Hinshaw I) have httle to no effect on 
asymmetry. 

Planetary Cut. If a planet is determined to be within ^cut = 1-5° of either radiometer's 
normal vector, the temperature map is not updated for either. This value is conservative, 
being over four times the FWHM of the V-band radiometer beam size (this helps mitigate 
the effect of non-Gaussian shoulders in the beam profile). We test the effects of making maps 
assuming 6'cut = 0.5°, and ^^cut = 0.0°. In both cases, asymmetry becomes more pronounced 
over the full range / [2,64]. If 6'cut = 0.5°, the effect upon significances is < 10%, while 
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Fig. 13.— Plots of c^giobai as a function of shifted dipole velocity (Ic = 1.87 km s for 
the multipole ranges displayed in each panel. The dot-dashed line in each panel represents 
Q^giobai = 0.05. This figure illustrates that changing the dipole velocity can have a marked 
effect upon the significance of asymmetry at the largest scales. However, since the dipole 
velocities used to generate this figure differ from the dipole velocity used to cahbrate the 
data, the uncertainty of each datum is unknown. 
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Fig. 14. — Uncorrected Kp2 cut-sky power (circles) and uncorrected Kp2+hemisphere cut- 
sky power (triangles and boxes for north and south respectively) for the dipole (top), 
quadrupole (middle), and octopole (bottom), as a function of change in the magnitude 
of the dipole velocity {la = 1.87 km s^^). The assumed NP is the galactic NP. This figure 
shows that changing the dipole velocity acts to increase the Kp2-|-hemisphere cut-sky power 
for the quadrupole, so that asymmetry approaches unity and becomes insignificant. It also 
demonstrates that for the Kp2 cut-sky, the effect of shifting the velocity is concentrated at 
the dipole, as expected. 
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Fig. 15. — Same as the upper left panel of Figure 10, with the magnitude of the dipole 
velocity increased by la (= 1.87 km s~^), 2a, and 3a (upper right, lower left, and lower right, 
respectively). The effects of increasing the dipole velocity are to make ctgiobai insignificant, to 
make the direction of maximum significance migrate toward the ecliptic plane, and to make 
the apparent great circle less visually striking. 
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if ^cut = 0.0°, CKgiobai — 0. The former result indicates that reasonably changing ^cut would 
have virtually no effect on asymmetry. 

Foreground Maps. We test the effect of foreground corrections by changing the relative 
normalizations of the three foreground maps (dust, Hq, and synchrotron; see Bennett II for 
details). We find that small shifts (~ 10%) have little to no effect upon asymmetry, which 
is not unexpected because of the markedly reduced level of foreground emission outside 
the Kp2 mask with respect to that in the galactic plane. (Because of the large number 
of normalization combinations tested, we do not display our results in Table 1.) To have 
a discernible effect, order-of-magnitude changes in the normalizations are required. This 
demonstrates that asymmetry is not simply related to uncertainties in relative normalization. 

5. Summary and Conclusions 

In this work, we analyze first-year WMAP data to determine the significance of asym- 
metry in summed power between opposite arbitrarily defined hemispheres. We perform this 
analysis on maps that we compute directly from the calibrated time-ordered data, using a 
map-making algorithm that we have developed to be similar to that used by the WMAP 
team. By creating our own maps, we are able to determine whether altering elements of the 
map-making algorithm is sufficient to affect the observed asymmetry (and perhaps to render 
it statistically insignificant). 

We follow Eriksen I and Hansen I by analyzing co- added V- and W-band foreground 
corrected maps, and like these groups we find that there is a significant difference in summed 
power between arbitrarily defined northern and southern hemispheres. We compute signifi- 
cance by simulating data from the running-index LCDM power spectrum, determining the 
most significant deviation from expectation in each dataset, fitting an extreme- value dis- 
tribution to these deviations, and computing the tail integral of the best-fit distribution. 
Over the multipole range I = [2,64], we estimate the global significance of asymmetry to 
be ctgiobai ~ 10""^, and we find that this value is not sensitive either to frequency (as deter- 
mined by examining co-added data in the Q, V, and W bands, and within each radiometer 
alone), or to the details of the power spectrum used as input to simulations. These results 
clearly indicate that the posited power asymmetry in WMAP data is real. We determine 
the smallest multipole ranges for which ctgiobai ^ 0.05; we find twelve such ranges, including 
I — [2,3] and [6,7], for which Ogiobai 0. Contour plots made for these ranges indicates 
an association between the direction of maximum asymmetry and the ecliptic plane (four 
of twelve directions within a; 4° and eight within a; 20°; Pbinomiai % 0.01), although one 
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must be mindful of the coarseness of our coordinate grid. At large scales, the contours of 
least significance (odir ~ 0.2-0.5) may be fit by eye with great circles inclined relative to the 
ecliptic plane, with inclinations of ~ 30° and 90° for / = [2,3] and [6,7] respectively. We find 
that the / = [2,3] great circle is insensitive to frequency and that it passes over preferred 
quadrupole and octopole axes derived by de Oliveira-Costa et al. and Schwarz et al., as well 
as the pointing direction of the CMB dipole. 

We test the robustness of our results by using examining the effect of map-making algo- 
rithmic assumptions upon asymmetry. We create and analyze new maps after changing the 
monopole temperature, the magnitude of the dipole velocity, the masks, the method of inter- 
polation used to determine the pointing direction of radiometer normal vectors as a function 
of time, the planetary cut, and the relative normalizations of the foreground maps. In vir- 
tually all cases, the alterations had little to no effect upon asymmetry either for the range 
I — [2,64] or at the largest scales. An exception to this is the magnitude of the dipole velocity: 
we find that increasing it has a marked effect on the direct (i.e., uncorrected) estimate of 
quadrupole power in the Kp2+hemisphere cut skies, increasing power in both the northern 
and southern hemispheres such that asymmetry approaches unity (and becomes insignifi- 
cant). Changing the magnitude of the dipole velocity is a frequency-independent change, 
consistent with observations that the significance of asymmetry is frequency-independent. 
If the true dipole velocity magnitude is ~ 1-3(7 larger than the current value, large-scale 
asymmetry may disappear. However, it will not influence asymmetry estimates at smaller 
scales, nor appreciably change full-sky power spectra. This is consistent with the conclusion 
of Hansen I that whatever causes asymmetry at large scales is not what causes it at smaller 
scales; however, they posit that galactic contamination is the cause of large-scale asymmetry. 

The effect of changing the magnitude of the dipole velocity, along with the observation of 
structures in our contour plots that follow great circles inclined relative to the ecliptic plane, 
suggests that the observed power asymmetry may be (at least partially) caused by the use 
of an incorrect dipole vector in combination with a systematic or foreground process that 
is associated with the ecliptic. The natural next steps are to determine the dipole velocity 
vector that minimizes asymmetry (with error bars), to make new maps based on that vector, 
and to examine the frequency dependence of any remaining significant asymmetry to attempt 
to differentiate between systematic or foreground causes. However, without the ability to 
recalibrate the first-year WMAP data, such an exercise may not yield accurate results. We 
thus would ask that both raw time-ordered data and calibration software be made available 
to the public in future data releases. 
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A. Map-Making Recipe 

In this Appendix, we lay out the details of our map-making algorithm, summarized in 

§2. 

A.l. Preliminaries 

We begin by assuming two values: the velocity of the Sun relative to the CMB rest frame 
in galactic coordinates, ("^Q^gab "^q gau "^©.gai) ~ (—26.26,-243.71,274.63) km s^^ (Bennett I), 
and the monopole temperature Tq = 2.725 K (Mather et al.). We input the planetary 
cut angle ^cut (default 1.5°), the method of quaternion interpolation (none, linear, or using 
Lagrange interpolating polynomials, with the latter the default), and the number of iterations 
(default 20). 

We input the mask from wmap_composite_mask_yrl_vl .f its. The masking column in 
this file is bit-coded; we examine bit 5 to set the processing foreground mask, which is the 
Kp8 mask with no edge smoothing. We do not apply the point source mask, as we expect 
point sources to have httle effect on results for Z < 64. 

We then read in a list of galactic coordinates for each HEALPix pixel, generated once 
using the IDL routine healpix_nested_vectors. We assume nside = 512 (or 3,145,728 
map pixels). 
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A.2. Mciking Map Estimates From Time-Ordered Data 

The making of maps involves four embedded loops: looping over iteration number, 
looping over the 366 time-ordered data (TOD) files/^ looping over the science frames within 
each TOD file, and looping over the data vector contained within each science frame. 

During each loop, we open and read each TOD file, as opposed to keeping their contents 
in memory after a first reading. (This is not strictly necessary; see §A.4 below.) Each TOD 
file contains five tables, of which we are interested in three: the Meta Data, Line-of-Sight 
(LOS), and Science Data tables. 

Each Meta Data table contains 1875 rows (recorded every 46.08 s, or 30 science frames). 
We extract from this table elements of the time column, the spacecraft position column 
(needed for determining planetary positions), the spacecraft velocity column, and the column 
containing the 33 x 4 quaternion matrices that encode the spacecraft orientation for each 
of science frame (the extra 3 elements provide seamless interpolation; see below). The 
spacecraft position and velocity are given in celestial coordinates; they are converted to 
galactic coordinates as detailed below. 

The LOS table contains the orientation of each radiometer horn in spacecraft coordi- 
nates, which we extract once. This information, combined with the information from the 
quaternion matrices, allows us to determine the normal vector for each radiometer horn in 
galactic coordinates, and thus the map pixel numbers associated with each radiometer as a 
function of time. We return to this below. 

Each Science Data table contains 56250 rows (30 x 1875), corresponding to science frames 
of length 1.536 s. Prom each row, we extract a data vector of length A^data (=15 for Q-band 
data, 20 for V-band data, and 30 for W-band data), and a status code. If this code is odd 
(i.e., if bit 1 is set), the data within the vector are considered bad and are not used. 

For each datum d, an interpolated quaternion is computed, using a Lagrange interpo- 
lating polynomial formulation provided by the WMAP team (G. Hinshaw & P. Butterworth, 
private communication). This formulation requires the quaternions for the previous, cur- 
rent, and following two science frames. We calculate the time At between that of the 
datum and the beginning of the current science frame, and convert this to a fractional offset 
/ = At/1.536 s. We then define four weights: 

w, = _/(/_i)(/_2)/6 
W2 = (/ + l)(/-l)(/-2)/2 



http : / /lambda . gsf c . nasa.gov/product/map/drl/map_tod . cf m 
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^3 = -(/ + l)/(/-2)/2 
W4 = (/+l)/(/-l)/6, 

and determine the four elements of the interpolated quaternion: 

4 

gint,i = YwjQj^i, (Al) 

i=i 

where j represents a science frame. The next step is to normalize g'int and to convert it to a 
transformation matrix (here we drop the "int" subscript): 

?1 - ?2 - ?3 + ?4 2(^1 g2 - ?3?4) 2(^153 + g2?4) \ 

2(?ig2 + g3?4) -g? + gi-g3 + ?4 2(g2g3-gig4) (A2) 
2(gi?3 - ?2g4) 2(52?3 + ?i?4) -g? - gi + gi + g| / 

This transformation matrix is in turn used to convert the radiometer horn normal vector Ugp 
from spacecraft to celestial coordinates: 

(^l,sp ^2,sp ^3,sp 

ni,sp ^2,sp ^3,sp I (A3) 

^l,sp ^2,sp ^3,sp 
3 

ni,cel = ^Nij^cel- (A4) 

We then convert to galactic coordinates using the transformation matrix 

(cos ^ cos $ — sin ^ sin $ cos sin ^ cos $ + cos ^ sin $ cos sin $ sin \ 
— cos \1' sin $ — sin ^ cos $ cos — sin ^' sin $ + cos ^ cos $ cos cos $ sin j , 
sin ^ sin — cos^sin0 cos0 / 

(A5) 

where * = 282.85948°, $ = 327.06808°, and = 62.871750°. The transformation is: 

^gal = Tcgncel- (A6) 

The mathematics of the subsequent transformation from figai to HEALPix pixel number is 
complex and will not be reproduced here; our software routine is based upon the HEALPix 
IDL routine vec2pix_nest. 

At this point, we know to which pixels the A and B radiometer horns are pointing. 
The next step is to determine the Doppler shift of the monopole relative to the spacecraft. 
The spacecraft velocity is recorded every 46.08 seconds, as noted above. Because the velocity 
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varies little during that amount of time, it is not strictly necessary to interpolate to determine 
velocities for each datum. (However, we have made the choice to linearly interpolate the 
velocities.) We convert the velocity associated with the datum from celestial to galactic 
coordinates: 

^''sp.gal Tcgt^sp^cel • ("^'^) 

Now we can remove the contribution of the Doppler-shifted monopole from the datum (here, 
we drop the "gal" subscript) 

/3 = ^{v.p + Vq) (A8) 



P-{fiA-fiB) + {P-nAy-{P-nBy , (A9) 



and update the temperature map and the counter recording how many times each pixel is 
observed: 

+ wt (AlO) 

Arp _ ^rj. ATcMB.t + (1 - a^im) ArcMB,pB,i ..^-.x 

^J-GMB,pA,i+l — ^J-CMB,pA,i+l + T— , \^^^) 

J- ~r 2^im 



and 



+ Wt (A12) 

.rp Arp I ... -ATcMB,f + (1 + a^im)ATcMB,pA,i /AioN 



where i + 1 is the current iteration (ATp^o is the input zero-tcmpcrature map). The counter 
is actually updated during only the first iteration as nobs,p does not change from iteration 
to iteration. As noted in §2, we assume the statistical weight Wf = 1.0 for all data, as the 
first-year data release does not contain information on the statistical weighting model used 
by the WMAP team. At the end of a given iteration, we divide by riobs.p to determine the 
weighted average. 

Equation (3) is skipped (i.e., the map is not updated) if pixel B lies within the processing 
mask, while equation (5) is skipped if pixel A is within the mask. Both equations are skipped 



^''What we denote ATcmb,« actually consists of contributions from the CMB and foregrounds. One can 
correct for the foregrounds during mapping by subtracting Tj^ _ T^'^ and (if desired) Tl'^0sp ■ ^a) - 
^B^(Ap ' "b) from ATcai.t, along with the contribution of the Doppler-shifted monopole. The latter factor 
assumes that the foreground emitters are stationary in the barycenter frame of reference. We find that the 
difference between performing foreground corrections during mapping and after mapping is negligible (< 1 
/xK in each sky map pixel). 
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if either pixel A or S is within an angle of ^cut degrees of any outer planet (except Pluto; 
see §A.3 below). 

Following Hinshaw I, we generally iterate 20 times; in Figure 16, we show that this 
number of iterations is sufficient for convergence. 

The map that is generated after 20 iterations is not an absolute map, but a differ- 
ence map with no precisely defined zero-point. The WMAP team, which like us uses a 
zero-temperature input map, determines the zero-point by fitting a function of the form 
T — Cicscfegai + Co to southern galactic hemisphere data below 6gai = —15° for each map 
(Bennett I). We apply a constant offset to our maps such that the mean pixel temperature 
matches that of the WMAP maps. 



A. 3. Planetary Cuts 

In theory, we should determine the positions of all outer planets in each science frame. 
However, position determination is a major computational bottleneck. Hence whenever 
we determine planetary positions, we estimate the length of time necessary for WMAP to 
rotate through the smallest angle between either radiometer normal vector and any outer 
planet, and we do not determine positions for that length of time. This scheme is unduly 
conservative, since we do not check to see if WMAP is rotating toward that planet, but it 
effectively eliminates the bottleneck. 

We compute planetary positions using the publicly available axBary software package^^. 
Wc use the package function dpleph, based on the cphcmcris JPLEPH.405. For each outer 
planet, wc (a) determine its position relative to the spacecraft in celestial coordinates for 
Julian date tjD, applying spacecraft positions read in from the Meta Data table; (b) update 
that position given the actual time of the science frame; (c) convert that position to galactic 
coordinates using the transformation matrix Teg; and (d) calculate ^a/b = acos(npianet ■ 
fiA/B,gai)- Note that we do not include the effect of light travel time in our algorithm, as 
the planets move by less than one sky pixel (~ 7') during that period. If 9x or is < ^cut 
(default 1.5°), we do not update the temperature map for either pixel. 



http : / /heasarc . gsf c . nasa . gov/listserv/heaf it s/msg00050 . html . 
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Fig. 16. — The average change in temperature in each map pixel from one iteration to the 
next, as a function of iteration, for all pixels (top) and for pixels outside the Kp8- without- 
edge-smoothing mask (bottom). The particular values shown here were generated for a 
40-iteration Wl radiometer map. The la error bars indicate the size of the distributions of 
AT, and are not estimates of the errors on the mean. These plots demonstrate the robustness 
of qualitative conclusions based on 20-iteration maps. 



A. 4. Computational Details 



Map-making is a computationally intensive task. Given current typical CPU speeds 
(~ GHz), making a single 20-iteration map takes ~ 1 CPU day. To reduce map-making 
time, we parallelized our C code using the Message Passing Interface (MPI) library. In 
this instance, parallelization is particularly easy: each processor is assigned to loop over a 
subset of TOD files, taking the estimated map from the previous iteration and updating it 
only for that subset. At the end of each iteration, the updated maps are broadcast to the 
root node, which sums them; the summed map is then broadcast back to the processing 
nodes for use during the next iteration. We run our parallelized code on the Teragrid Linux 
Cluster (www. teragrid. org), which links massively parallel super computing clusters at nine 
sites, including NCSA, where we do the bulk of our map-making and where a typical 24- 
processor, 20-iteration run takes ~ 25 CPU minutes. We note that although we read in the 
data sequentially during each iteration, because each processing node has its own memory 
we could in theory keep the time-ordered data in memory after being read once. The total 
amount of, e.g., W-band data that would need to be kept in memory would be reduced from 
366 days x 56,250 data rows x 30 data = 6.2 x 10^ (r^ 2.5 Gb of memory) to, say, 2.6 x 
10^ (Ri 100 Mb of memory) for 24 processors, an amount that each node (containing two 
processors) could easily handle. 



http : / /www-imix . mcs . anl . gov/mpi 
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